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in Abstract 

We compute El transitions and electric radii in the Beryllium- 11 nucleus using an effective 
'"^ field theory that exploits the separation of scales in this halo system. We fix the leading-order 

4^ parameters of the EFT from measured data on the 1/2"*" and 1/2^ levels in ^^Be and the B(E1) 

strength for the transition between them. We then obtain predictions for the B(E1) strength for 
^ Coulomb dissociation of the ^^Be nucleus to the continuum. We also compute the charge radii of 

I ^ I the 1/2"'" and l/2~ states. Agreement with experiment within the expected accuracy of a leading- 

order computation in this EFT is obtained. We also discuss how next-to-leading-order (NLO) 
^ corrections involving both s-wave and p-wave ^*^Be-neutron interactions affect our results, and 

C display the NLO predictions for quantities which are free of additional short-distance operators at 

00 

this order. Information on neutron-^'^Be scattering in the relevant channels is inferred. 
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I. INTRODUCTION 



The first excitation of tlie BerylIium-10 nucleus is 3.4 MeV above tlie ground state, 
and tliat ground state lias spin and parity quantum numbers = O"*". Meanwhile, the 
Beryllium-11 nucleus has a 1/2+ state whose neutron separation energy is 500 keV, and a 
l/2~ state whose neutron separation energy is 180 keV [1]. The shallowness of these two 
states of ^^Be compared to the bound states of ^°Be suggests that they have significant 
components in which a loosely-bound neutron orbits a ^°Be core. In this "one-neutron halo" 
picture the 1/2"*" is predominantly an s-wave bound state, while the l/2~ is predominantly 
a relative p-wave between the neutron and the core. In this paper, we discuss efforts to 
use effective field theory (EFT) to systematically implement such a halo picture of the ^^Be 
nucleus. 

This halo viewpoint is reinforced by the fact that the scattering volume of n-^°Be scat- 
tering in the / = 1, J = 1/2 channel has been determined to be [2] 

ai = (457 ± 67) fml (1) 

The corresponding length scale of order 8 fm is large compared to the natural length-scale 
of core-neutron interactions, which is ~ 2-3 fm. 

The datum ([T]), together with the information on the bound-state energies in the ^°Be 
and ^^Be systems, helps us to estimate the expansion parameter in our Halo EFT. This is the 
binding energy of the halo nucleus, as compared to the energy required to excite the core, 
i.e. i?io/-Bhi ~ 1/6. Converting this to an estimate of the different distance scales involved, 
we infer that a majority of the probability density of ^^Be occupies a region outside the ^°Be 
core: -Rcore/-Rhaio ~ 0.4, which is consistent with the ratio implied by the numbers in the 
previous paragraph. This ratio of distance scales is the formal expansion parameter for the 
EFT, and since it is not particularly small, leading-order calculations are only a first step. 
We therefore present calculations up to next-to-leading order for several quantities, in order 
to confirm that the series is converging as expected. 

In particular we apply this EFT to electromagnetic reactions in the ^^Be system. The 
B (El) (1/2+ — !■ l/2~) transition has recently been measured to be 

B(El)(l/2+ ^ 1/2-) = 0.105(12) fm^ (2) 

using intermediate-energy Coulomb excitation This is consistent with the older number 

B(El)(l/2+ ^ 1/2^) = 0.116(12) fm^ (3) 

from lifetime measurements [1] . There are also two recent data sets on the Coulomb-induced 
breakup of the ^^Be nucleus [5l |6] (see also Ref. [7J). Both experiments extracted the 
excitation function dB{El)/dE as a function of the energy of the outgoing neutron E. 
For low neutron energies this excitation function is affected by the final-state interaction in 
the p- waves, and can be predicted in the halo picture [2]. Ref. [5] also extracted a neutron 
radius for the ground state of ^^Be from their data: 

(r2)i/2 = 5 7(4) fjn. (4) 
This is consistent with the recent atomic-physics measurement of the ^^Be charge radius [8]: 

{rDlil = 2.463(16) fm. (5) 
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All of these measurements can be addressed within the Halo EFT we will use here. In this 
theory the s- and p-wave states of the Beryllium-11 nucleus are generated by core-neutron 
contact interactions. The theory does not get the interior part of the nuclear wave function 
correct, but, by construction, it reproduces the correct asymptotics of the wave functions of 
these states: 



for the 500 keV and 180 keV states, respectively. The quantities 70 and 71 are determined 
by the neutron separation energies of the states in question. At leading order (LO) in the 
expansion the Asymptotic Normalization Coefficients (ANCs) Aq and Ai are fixed. (In the 
case of the p-wave this is related to the theorem discussed in Ref. [9].) However, at next- 
to-leading order Aq and Ai become, in essence, free parameters of the theory, and must 
themselves be extracted from data. 

Halo EFT is well-suited for this task. It is not intended to compete with ab initio calcula- 
tions of this halo nucleus (see, e.g. (TUl E]) or of ^''Be-n scattering [12], or with microscopic 
descriptions of the ^^Be El strength (see, e.g. Instead, Halo EFT is complementary to 
such approaches, since it takes Aq, Ai, 70, and 71 as input, rather than seeking to predict 
them via a detailed description of the A = 11 system. The EFT's goal is to ensure that the 
long-distance properties of the halo are correctly taken care of, and then to elucidate the re- 
lationships between different observables in the ^°Be-n system that result. Such correlations 
fiow directly from the existence of the shallow 1/2"*" and l/2~ bound states in this system. 
In particular, below we will show how Aq and Ai are correlated with neutron-^^Be scattering 
observables, as well as with B(El)(l/2+ — )■ l/2~) and the Coulomb dissociation data. This, 
in turn, demonstrates how — and with what accuracy — the ANCs can be inferred from these 
experimental quantities. 

A preliminary version of these findings appeared in Ref. [I3]. The presentation here 
corrects and expands upon that earlier work. A related study of radiative neutron capture on 
Lithium-7 was recently carried out in [T3] . The mechanism for the cancellation of divergences 
in the s-to-p transition in this reaction and in the Beryllium-11 system is the same. 

II. HALO EFT FOR BERYLLIUM-11 

We use the "Halo EFT" developed in Refs. [151 [IS] to calculate the properties of the 
Beryllium-11 nucleus. The degrees of freedom in our Halo EFT treatment are the ^'^Be core 
and the neutron. The EFT expansion in this case is an expansion in powers of w/Shigh- 
Here -Bhigh is, e.g. the excitation energy of states in ^°Be, and so is of order a few MeV, and 
u is the energy of the photon exciting the electromagnetic transition of interest. 

A. Lagrangian: strong sector 

In our calculation, the ^°Be core and the neutron are represented by the bosonic field c 
and a spinor field n, respectively. We include the strong s-wave and p-wave interactions that 
lead to the shallow bound states in the ^^Be system through the incorporation of additional 



Uo{r) 



Aq exp(-7or) 




(6) 
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fields, which encode physics of these states. Therefore the 1/2"'" state is constructed as a 
spinor field, a^. In contrast, the field representing the l/2~ state is a pseudo-spinor, i.e. it 
is parity odd. We denote it as vr^. Its behavior under parity restricts the types of couplings 
which can appear in the Lagrangian. In particular, only combinations of nucleon and core 
fields with an odd number of derivatives can couple to vr. Moreover, when constructing these 
operators we must decompose them into their irreducible representations under the rotation 
group. The portion that couples to vr^ is then just the J = 1/2 part. For example, if we 

construct nj3{i Vj)c where Vj= (V — V)j, this combination has both J = 1/2 and J = 3/2 
parts. We project out the J = 1/2 part by defining: 



[n{i V)c]i 



/9i 



Up 



•f-> 



(7) 



where {jimij2m2\{jij2)JM) is the Clebsch-Gordan coefficient to couple ji and j2 to J. 

With the masses of the neutron and the core are denoted by m and M, respectively, and 
Mnc = {M + m) is the total mass of the n-^'^Be system, we then have: 

V2 



c 



+CTI 



tdt + 



2M 



-go [c"^-"^ 
91 M 



c + n 
2M„, 

rt. 



Ao 



Vi ( idt + 



2M„ 
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&nlas + cringe] — 

i V {nc 



m 



2 M 



7r 



n{i V)c ^ 



+ 



+ + 



Ai 



7r< 



TT, 



(8) 



where we adopt the convention that repeated spin indices are summed. Note that the last 
line of Eq. ([s]) represents an additional p-wave interaction necessary to maintain Galilean 
invariance. It is required because the n and c fields have different masses. The dots represent 

higher-order interactions not considered here. One such interaction involves the J = 3/2 

•f-> 

part of the [n{i V)c] operator: 



3/2 



C(3/2) 
C(3/2) 



asijs2k 



1 



1 



-a 



-S2lk 



-a 



k)c) 



(9) 



As we shall discuss in the next section, this interaction is assumed to be natural in 
our power counting, in contrast to the interactions mediated by vr and a fields, which are 
enhanced. 



B. s-wave Be- neutron interactions 

In order to treat the shallow s-wave state in the ^°Be-neutron system we adopt the 
counting that has been successfully developed to treat shallow s-wave states in the nucleon- 
nucleon system [T7H2T] . In leading order, the a field is static and its bare propagator is simply 
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FIG. 1: Diagrammatic representation of the Dyson equation for the dressed a propagator repre- 
senting the s-wave ^^Be-n bound state in the theory. Here and below the dashed hne indicates the 
field for the ^"Be core, and the thin solid line is the neutron. The thick grey line is the bare a 
propagator, and the thick black line is the dressed a propagator. 



FIG. 2: The Feynman diagram relating the dressed a propagator to the s-wave neutron-core 
scattering amplitude. 

I/Aq. Dividing out the leading term, the one-loop correction to the bare propagator is g"^/ Aq 
times an nc bubble which has a typical size of order l/i?haio ^- Since Qq/Aq ~ -Rhaio, this 
correction is of order one. Consequently, the nc loops must be resummed when computing 
the full a propagator. 

This can be achieved through the Dyson equation shown in Fig. [T| which leads to: 



Ao + Vo[po - pV(2M„,) + ie] - S,(p) 



(10) 



with So-(p) the one-loop self-energy for the a field. 
This one-loop self-energy is calculated as: 



glrriR 
27r 



2M.. 



+ /U 



when computed in power-law divergence subtraction (PDS) with a scale n 
have introduced the reduced mass of the neutron-core system: 



(11) 



Here we 



mM 



rriR 



m + M 

and the limit e — 0^ in the end is understood. 



(12) 



^ In a suitable regularization scheme, e.g. power-law divergence subtraction |18l I19| . this is true for both 
the real and imaginary parts of the loops. 
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Substituting Eq. (11) into Eq. (10), we can set the parameters 5^0 and Aq by computing 



the s-wave neutron-core scattering amphtude in the theory defined by Eq. ([8]) (see Fig. [2]): 

t,{E)=glD,{E,Q), (13) 

in the two-body center-of-mass frame with E = A;^/(2m/j). This is then matched to the 
efi^ective-range expansion in this channel: 

2n 1 



to{E) 



producing 



rriR 1/ao - ^rok^ + ik 
27r7o 1 1 



^r9o 1 - rolo po 



2M„ 



+ Ra{p) 



(14) 



(15) 
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where Ra{p) is regular at the pole po — = —Bq. In Eq. (15), the position of the pole 
is determined by the binding energy Bq = '~^l/{2mji), and 70 is the positive root of the 
equation: 

1 1 



+ o'^oTo - 7o = , 



(16) 



The wave-function renormalization for the dressed a propagator can be read off as the residue 

ry 27r7o , 



of the pole in Eq. (15): 



(17) 



This result is valid to NLO in the expansion in Rcore/Rhaio ~ ^olo- It yields a wave function 
(|6| with 



An 



270 



(18) 



C. p-wave ^''Be- neutron interactions 



We proceed similarly for the p-wave state. The propagator, -D^(j9)ss/ = Dt,{p)Sss', for this 
state obeys the Dyson equation depicted in Fig. [3) Rather than computing the self-energy 

of the TT field directly it is easier to compute the self-energy for a field 7^/3(2 Vj)c and then 
couple the neutron spin and the relative momentum in the appropriate way to project out 
the J = 1/2 piece of the result. Hence we now consider the one-loop self-energy, S7r(p)ij,/3a 
for such a p-wave field. We first observe: 

S^(p)ij>/3 = SijSapT.{p) . (19) 

The scalar function: 



E(p) 



6tt 



'^rriR Po - 



2M„ 



2m/j Po - 



2M„, 



+ te ] + fi 



(20) 



where the PDS scheme has been employed and momentum traces have been performed in 
three dimensions. From this we can construct a self-energy for transitions from the vr-field 
state s to the 7r-field state s': 



sip) 



S(p) 



(21) 
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FIG. 3: Diagrammatic representation of the Dyson equation for the dressed vr propagator rep- 
resenting the p-wave ^''Be-n bound state in the theory. Once again, the dashed hne indicates the 
field for the ^"^Be core, and the thin sohd hne is the neutron. The thin double line is the bare tt 
propagator, and the thick double line is the dressed vr propagator. 



since S is independent of /3 and j we can use completeness of the Clebsch-Gordon coefficients 
to show that S is diagonal in s' and s, i.e. S^3;j,(p) = Sg's'^ip)- It follows that D^^lp) takes 
the form: ^ 

^ Ai + r/ibo-pV(2M„c)]-S(p) • ^^^^ 

We note that since the self-energy loop is cubically divergent both parameters, Ai and 
Qi are mandatory for renormalization at LO This time we are interested in the p-wave 
core-neutron scattering amplitude in the center-of mass frame: 



ti(p',p;E) = gfp' ■pD^{E,0) 



P - P 



niji 1/ai — ^rik^ + ik^ 



(23) 



with k = ^/2/nv[lE = \p'\ = |p| for on-shell scattering. Consequently we obtain: 



Qtt 



1 



+ regular . 



(24) 



Here 71 = y/2mfiBi is the solution of 



11.0 
- + + 7? = 



(25) 



where ai is the scattering volume, and ri the p-wave "effective range", which, in fact, has 
dimensions of 1/length. Both parameters are required to leading order in the Halo EFT. The 
wave-function renormalization for the dressed tt propagator can be read off as the residue of 



the pole in Eq. (24): 



1 



"^Rfl'l ^1 + 371 



(26) 



The propagator ( 23 ) has three poles corresponding to th e zero es of Eq. ( 25 ) . Using the 
NLO parameter values for the ^°Be-n system obtained in Sec. H F below, we find two bound- 
state poles corresponding to typical momenta 71 ~ l/i?haio and 71 ~ l/i?core- The first is 
that which we identified with the ^^Be excited state in the previous paragraph. The second 
is a spurious bound state, which is outside the domain of validity of Halo EFT, and is not 
physically realized in the ^^Be system. The third solution of Eq. (25) represents a virtual 
state with a typical momentum ~ l/i?haio- 

The power counting for the propagator Dt^{p) that we adopt here is that of Ref. [T6] . 
We take ri ~ l/i?core- The propagator then has a pole at 71 ~ l/-Rhaio5 which occurs 



7 



"kinematically" when rifc^ ~ 1/ai. ai then must obviously be of order R\^Yo^core 



for such a 



"kinematic" pole to occur. We note that for k ~ l/-Rhaio the unitarity piece of the propagator 
has a size Thus, away from the pole, the dominant contribution to the vr propagator 

now comes from the bare part (after appropriate renormalization) and so: 

Gtt P • P 



ti(p',p;E) 



(27) 



where we have also used Eq. (25) to re-express 1/ai in terms of 71 and then dropped the 7^ 



term relative to the (larger) r^Tf piece. 



The result (27) can be easily re-expressed as: 



P P 



myri E + Bi 



(28) 



The amplitude (28) has only a pole at E 

to the 1/2- state of ^^Be 



— Bi = —jf/{2mji) — a pole that corresponds 
(This pole actually occurs on both sheets of the complex E- 



plane, since (27) exhibits poles in both the lower and upper half of the complex fc-plane.) 



In contrast, the spurious deep pole from Eq. (23) has disappeared from the expressions (27) 
and (28). Our power counting therefore reproduces the spectrum of the ^^Be system. The 



p-state wave-function ^ is then obtained, with 



27? 
-ri' 



The p-wave phase shifts in this theory are given by: 

1 



k"^ cot 5i 



(29) 



(30) 



if no expansions are made. However, for ri ~ l/i?core we can, once again, drop 75* to leading 
order in a power counting in Rcore/ Rhaio, with the result that 



cot 61 



ri 
2 



k k^ 



+ 



Rc 



R 



halo 



(31) 



Since ri ^ fc, 71 we have cot 5i large, which implies that 5i is approximately zero. Indeed 



5i 



k^ 



+ 



Rr 



R 



halo 



(32) 



ri k'^ + 7^ 

for all A; ~ 7i ~ l/i?haio- Small phase shifts imply small unitarity corrections, which is why 
the imaginary part of t^^ can be treated perturbatively in this regime. 

The only exception to this occurs if we consider \E — Bi\ ~ In that case we are 

in Eq. (27) cancel, or come close to 
k^ and 



close to the pole and the two terms ~ (-Rc 
doing so. It then becomes necessary to resum the pieces 



"-core-Rhalo) 



'jf which were dropped 

in order to obtain Eq. (28) [151 [22] • In particular, if _Bi < (i.e. the pole is at positive 
energy) then these corrections shift the pole off the real axis and mean that it represents a 
resonance. The propagator ( [23] ) thus describes a p-wave resonance if the pole is at positive 
energy and the width of the resonance at energy Er is of order k^/ tiEr [15l UHl |23]. This is 
thus a narrow resonance if ri ~ l/i?core- The case of p-wave resonances will not be discussed 
further here, since we will restrict ourselves to Bi > 0, as is relevant for ^^Be. This is the 
case of a shallow p-wave bound state. 
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D. Rescaled fields and naive dimensional analysis 



At this point it is useful to rescale the piece of the Lagrangian which encodes the enhanced 
s- and p-wave interactions. We rewrite L in terms of fields with non-canonical dimensions, 
which absorb factors of g^^ gi, m, and M pij. We define: 



= (^sgomR, TTs = -KsQimR . (33) 

The scaling of the absorbed factors goruR and giniR can be obtained by recalling the matching 
between (13) and (14) for s- waves, and in Eq. (23) for p-waves. At leading order, this yields: 



(34) 



In our counting, we have tq ~ 1/ri ~ -Rcore, which then determines how go and gi scale with 
i?core- Note also that, since ro > and ri < for n-^°Be interactions, we have rjo = — 1, 
r^i = 1 in this system. 

We analyze Mnc^ because this operator is dimension 5, and has no powers of the mass 
in it anymore [25] . The expression for the product of the total mass and the pertinent piece 
of the Lagrangian then becomes, in terms of these fields: 



^ at [r/o (2M„,9t + V") + M„,Ao] a. 



1 



fcl [r/i {iM^A + V) + M„,Ai] 



2mR 



^1 r + ^ + 

n{i V)c + c\i V)rv 



7r< 



(35) 



where the pieces of C that restore Galilean invariance have been suppressed. In terms of 
these new fields all the coefficients — even those in the "enhanced" interactions which generate 
shallow bound states (and resonances) are natural. The shallowness of these states is now 
encoded in the fact that the fields associated with them have non-canonical dimensions: 
[a] = 2, [vf] = 1, and non-canonical wave-function normalization — even at tree level. 



E. Lagrangian: electromagnetic sector 

Photons are then included in the Lagrangian via minimal substitution: 

d^^D^ = d^ + leQA^. (36) 

The charge operator Q takes different values, depending on whether it is acting on a c field 
or an n field. Qn = for the neutron, and below we denote the eigenvalue of the operator 
Q for the c field as Q c = QcC. Qc = 4 in the case of interest here, where the core is 
Beryllium- 10. 

Here our focus is on electric properties (and form factors), and the dominant pieces of 
the electric response can be derived by looking at how the Lagrangian ([8| is affected by the 



substitution (36). But, at higher orders in the computation of these properties, operators 



involving the electric field E and the fields c, n, a, and n which are gauge invariant by 
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themselves contribute to observables. Possible one- and two-derivative operators with one 
power of the photon field are: 



EM 



ll'j 

-LSvr|(V'Ao-9o(V-A))7r; 



3/21' 



ll'j 



2 2 



2 2 



(37) 



Note that if magnetic properties are to be considered we would also include operators in- 
volving diAj — djAi and the neutron, core, and bound-state fields. 

The electric interactions in Eq. (37) are gauge invariant by themselves, and so we must 



determine the order at which they occur. To do this we rewrite the Lagrangian (37) in terms 



of the rescaled fields (33). In terms of these fields we assume scaling by naive dimensional 
analysis with respect to the scale -Rcore of the operators that appear in MncC We then 
obtain the following scaling of the coefficients written above: 



'El ~ 



L 



(3/2) 
El 



r?3 7(o')„2_2 
-"'core'COyo''''iJ ' 

p 7(1/2)^ ^ ^2 
J^coret-El yoyi'''_R : 

p 7W^2 2 
J^core''COyi"''R i 

p4 /(3/2) 



(38) 
(39) 
(40) 
(41) 



where the parameters • are all of order one. 

Below we show that the leading effects in the El(l/2+ — t- l/2~) matrix element have 

Including the proper wave-function renormalization 



parametric dependence R 



halo 



-Rv, 



factors, the operator ~ L^U"^^ yields an effect ~ -Rcore a/ p 

-f^halo 



and so occurs already at NLO 



in that quantity. Similarly the leading effects in the charge-radius-squared of the 1/2 state 
in i^Be are ~ Rha\o/'ri ~ -Rhaio-Rcore- The operator proportional to L^Tg above produces effects 
in this charge radius of order i?core/^i ~ -Rcore) ^^^1 so affects the prediction for the p-wave 
radius at next-to-leading order. Thus if we desire NLO accuracy for quantities involving the 
shallow p-wave bound state there are two parameters in the Halo EFT description of the 
electric structure of the Beryllium-11 which cannot be fixed from ^'^Be-neutron scattering 
information alone. 



F. Fixing parameters 

Using the values Bq = 500 keV, Bi = 180 keV from Ref. pj, we infer 70 = 0.15 fm"^, 
and 7i = 0.09 fm"^, which are both of the expected size l/i?haio- From the power counting 



discussed in Sec. II C, we have at leading order: 



n = -—. (42) 

71^1 

It follows that if we adopt the value extracted in Ref. p] from experimental data, Eq. ([T]), 
we have ri = (—0.54 ± 0.08) fm"^. This number should, however, be taken as indicative 
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rather than definitive, since in the end we will fit ri to data at both LO and NLO, and 
deduce corresponding values for ai from our results. But already we see that the order of 
magnitude estimate implied by our counting ri ~ l/i?core is borne out by the numbers. 



At NLO, Eq. (42) is corrected to: 



ri = -— -271, (43) 

which, if we again use Eq. ([T| to get an idea of the effect, alters ri to (—0.72 ± 0.08) fm^-^. 
Such a ~ 30% correction is in line with the anticipated expansion parameter of Halo EFT 
in the ^^Be system. 

In the s-waves the situation is more straightforward: there we count ao ~ -Rhaio, and 



tq ~ -Rcorc- In consequence we can set tq = at LO, and obtain from Eq. (16) 



At leading order, in s-waves we have: 



70 = -■ (44) 

ao 



= ^ (45) 



and all other pertinent results can be obtained by taking the ro — )■ limit of the formulae 
in Sec. HTBI 

Thus, the parameters in Halo EFT for Beryllium-11 bound states at LO are ri, 70 (or 
equivalently gq), and 71. At NLO these are to be supplemented by r^, and the electromag- 
netic contact interactions for the 1/2+ — )■ 1/2^ El transition and l/2~-state radius. 



III. RESULTS FOR BOUND-STATE OBSERVABLES 



Using Eq. (|8j) plus minimal substitution (36), we obtain a Lagrangian that describes 
interactions amongst the core, the neutron, the ground and excited states of the ^^Be nucleus, 
and photons. In this section we use this Lagrangian to compute the form factor of the 1/2"'" 
and l/2~ states and the El transition from the 1/2+ to the 1/2^ state. 



A. s-wave form factor 



The s-wave form factor is computed by calculating the contribution to the irreducible 
vertex for A^aa interactions shown in Fig. |4} This is the only diagram it is necessary 
to consider at leading order. After the application of wave-function renormalization, the 
irreducible vertex for the Aq photon coupling to the a state is equal to —leQ cG E{\ci\) , where 
q is the three-momentum of the virtual photon. (Such an interpretation is valid provided 
the computation is carried out in the Breit frame, where the four-momentum of the virtual 
photon q = (0,q).) A straightforward calculation yields: 

G;^^(|q|) = |^-ctan({M) ' ^^^^ 
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FIG. 4: The LO contribution to the irreducible vertex for an Aq photon to couple to the field 
representing the ^"^Be-neutron s-wave bound state. Note that there is no diagram for the photon 
to couple to the neutron as this order, since Qn = 0. 



with / = m/Mnc = rnn/M. Note that Ge{0) = 1, as it should. For the deuteron, we have 
/ = 1/2, and this reduces to the LO result of Ref. [26] . 

The electric radius of the s-wave state can be extracted according to: 

4^^(|q|) = l-^(r|)^V + ..., (47) 



and an expansion of Eq. (46) in powers of |q| then yields 

{rir^ = {!-2- (48) 



Eq. (48) gives the electric radius of the ^^Be ground state relative to the electric radius of 
^°Be. In order to compare with the experimental radius, we therefore have to add our result 
and the charge radius of ^°Be in quadrature: 

(r|)nBe= (4)iOBe + f^. (49) 

This relation can be derived by writing the charge distribution of ^^Be as a convolution of 
the charge distribution of ^°Be with that of the ^°Be-n halo system. Using the convolution 
theorem for the Fourier transform, one finds that the total rms radius squared can be written 
as the sum of the squared radii for ^°Be and for the ^'^Be-n halo system. The latter effect 
can be calculated in the Halo EFT. Inserting the value 70 = 0.15 fm~^ obtained in the 
previous section and using the experimental result for the ^°Be charge radius [8], we find 
('"l;)iiBc ~ ^-40 from this leading-order HEFT computation. This is 2-3% smaller than 
the atomic physics measurement which yields {r\)\{'^^ = 2.463(16) fm [8]. In fact, comparing 
our result for (r^)iiBe — (?'c)i°Be (0.19 fm^) with the experimental result for this quantity 
(0.51(17) fm^), the agreement looks poor. But, this difference is actually consistent with the 
nominal 40% size of NLO effects when Halo EFT is applied to this system. 

At NLO a careful treatment of current conservation, which includes an operator associated 
with gauging the term ~ a'^doa in Eq. ([s]), still yields Ge(0) = 1, but also produces an 
increased charge radius, as long as tq > 0, cf. Ref. 



4)..Be=(r|).0Be + ^^^3{r^. (50) 
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Therefore NLO corrections improve the agreement with experiment. The precise size of the 
increase is fixed once the s-wave effective range tq is determined, as we shall do in Sec. fV] 
below. 

We can also obtain from our leading-order calculation a number for the mean-square of 
the relative core-neutron co-ordinate, (r^), i.e.: 

To convert this to a neutron radius, we must insert the conversion factor mc/{mc + rrin)- 
When this is done we find a LO neutron radius for the ^^Be ground state of: 

{riy/^ = 4.3 fm. (52) 

The neutron radius can be measured with a probe that couples only to the neutrons. To 
a very good approximation, the weak gauge boson Z° constitutes such a probe. Thus 
one could, in principle, measure the rms neutron radius using parity-violating electron 
scattering — c.f. the PREX experiment for the case of Lead-208 [2H]- However, a mea- 
surement of parity-violating electron scattering from ^^Be is certainly beyond present-day 
experimental capabilities. 



(51) 



B. p-wave form factor 



In this subsection, we calculate the charge form factor of the l/2~ excited state in ^^Be. 
NLO corrections might be expected to be smaller there since its binding energy, and so its 
typical momentum, is lower. However, as we shall see, a counterterm enters already at NLO 
in this observable. 

The p-wave form factor is computed by calculating the contribution to the irreducible 
vertex for Aotttt interactions shown in Fig. |5| There are two diagrams at LO. The first 
diagram is analogous to that for the s-wave state while the second diagram represents a 
direct coupling of the photon to the vr field. The latter contribution is leading order for the 
p-wave state since the effective range ri is leading order for this state. 

As with the self-energy of the vr-field, it is easier to compute the irreducible bubble for 

an Aq photon coupling to the p-wave field n(s{i Vj)c. In this way, we find that, after 
the application of wave-function renormalization, the irreducible vertex for the Aq photon 
coupling to the l/2~ state in the Breit frame can be written as: 



«(p')|JV.(p)) = -^eQ,5^(^al. (^l) (^Ij 

al3ij ^ V / / V 



1 \ 1 

-1 -s 
2/2 



S. 



Gi?'(|q|)*., + 



2M2 



(53) 



where q = p' — p is the three-momentum of the virtual photon. Here we have expressed the 
form factor in terms of the charge and quadrupole form factors of a vector field. Choosing 
q = |q|;z, and exploiting the fact that the neutron spin is unaffected by the charge operators 
that can occur up to the order we consider here, a brief calculation shows 



(7r;(p')|JV.(p)) = -^eQA'.G(,-H|q|). 



(54) 
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(a) (b) 

FIG. 5: The LO diagrams contributing to the irreducible vertex for an Aq photon to couple to 
the TT field representing the ^"^Be-neutron p-wave bound state. Diagram (a) is analog to the s-wave 
form factor while diagram (b) arises since the effective range ri enters at leading order for the 
p-wave state. 



The quadrupole form factor is thus unobservable in the l/2~ state. It could be observed in 
a 3/2~ state. The lowest 3/2~ state in ^^Be is 2.69 MeV above the ground state. However, 
this state is a n-^°Be scattering resonance that corresponds to typical momenta of order 
l/i?core and thus is outside the range of applicability of the Halo EFT. Since the situation 
could be different in other one-neutron halo nuclei, we quote the result for the quadrupole 
form factor for completeness: 



G'n\\ci\) 



2Ml 



n + 371 4|q|3/ 



2|q|/7i + (|qrr-47r)arctan 



m 

271 



(55) 



In the case of ^^Be, only the charge form factor is observable. A straightforward calcula- 
tion yields: 



G^?(|q|) 



1 



ri + 371 



ri 



|q|/ 



2|q|/7i + (|qrr + 27i)arctan 



M 

271 



(56) 



where again / = m/Mnc = itlr/M. For a strict LO result ri + 371 should be replaced by ri 
in these expressions. 

We have G^^\o) = 1, as required by charge conservation. The electric radius of the 
p-wave state relative to the ^°Be ground state can be extracted according to Eq. (47), and 
we obtain the electric radius for the case ri 
of |q|: 

5/^ 



2\W 



7i from an expansion of Eq. ( 56 ) in powers 

(57) 



271(371 + '^i) ' 

At LO in the situation of interest here |ri| ^ 71, we should replace this by: 



2\W 



5P 
27iri 



(58) 
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"'Core • 



This is, as promised above, ~ i?haio-Rc 

It is interesting to examine why the p-wave radius is more sensitive to short-distance 
physics than the s-wave one. From the co-ordinate space point of view we find that at LO 
the radius of the p-wave state must be calculated as: 

= Alf r drr^ (l + — y . (59) 

Jo V lif J 



Inserting Ai from Eq. (29) yields Eq. (58). Since the integral is finite, we can compute the 



contribution to it from values of r <^ 1 /7i : 

= Alf dr f 1 + — V e"^-- , (60) 
which produces, for -Rcore ^ I/ti, a short-distance fraction of the total result of: 

~ -rtcoreTl ■ (D-l-j 



The parametric dependence of the short- distance contribution on -Rcore is in accord with 
the result obtained in the previous section for the corresponding counterterm using naive 
dimensional analysis in our rescaled Lagrangian. 

It might seem counterintuitive that there is a short-distance contribution to {r'^)^'^^ al- 
ready at NLO — especially when the corresponding effect does not occur in (r l)^'^) until N^LO 



(see Eq. (38) and Ref. fSBj). The physics of this is, however, quite straightforward. It is 
associated with the propensity of the p-wave state's probability distribution to be drawn in 
to shorter distances than the s-wave one, as it gets caught between the attractive potential 
that produces the excited state of ^^Be and the centrifugal barrier. Observables associated 
with a shallow p-wave bound state will, therefore, generically exhibit counterterms at lower 
order than those of their s-wave counterparts. 



Numerical evaluation of the LO expression (58) leads to the prediction for the charge 
radius of the ^^Be p-wave state relative to the ^'^Be ground state {r'l)^'^^ = 0.36 fm^ at LO. 
(Here we have used as input the value ri = —0.66 fm, see the next section.) A trivial way 
to estimate the size of NLO corrections is to use the result for Ai in the limit |ri| ~ 71: 



A^ = \l %- (62) 

ri + 371 



instead of Eq. (29). This yields a 20% correction at NLO in agreement with the expectation 
from the power counting, although it must be remembered that the entire NLO result in 
the situation where |ri| ^ 71 includes the contributions of the operator ~ TrJV^AoVTs. Thus 
the only prediction we can offer here is a leading-order one. Using again the experimental 
result for the ^°Be charge radius [8], we predict the electric radius of the l/2~ state as: 

(4)i/L = (2.43 ± 0.1) fm (63) 

where our error bar comes from the above estimate of NLO effects. To our knowledge 
there is, as yet, no experimental determination of the charge radius of this state. Note that 
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scale with /, as 



in Eq. (63) we are assuming that the short-distance effects in 

suggestedby the renormahzation-group argument summarized in Eq. (60). Short- distance 
effects which do not involve a factor of /, e.g., modification of the proton distribution of 
^°Be in the ^^Be excited state due to a non-recoil effect, could have an impact on our final 
result that is larger than 0.1 fm. 



C. El transition: 1/2+ 1/2 state 



Now we consider the El transition from the 1/2"*" state to the l/2~ state. The irreducible 
vertex for this transition is depicted in Fig. |6} We compute the transition for a photon of 
arbitrary four momentum k = {u, k), and the sum of diagrams yields — ir^/j,^ where s' (s) is 
the spin projection of the 1/2^ (1/2^) state and /x is the polarization index of the photon. 




FIG. 6: The two diagrams needed for the irreducible vertex that governs the s-to-p-state transition, 
Tjfj_ in Halo EFT at leading order. 



We first observe that both the diagrams depicted in Fig. |6] are divergent, but that the 
divergences cancel, as they should since gauge invariance precludes us from writing down 
any contact interaction that contributes to this observable at leading order. Current con- 
servation at LO can be explicitly checked, and we find that, as long as both diagrams are 
considered [T5] : 

k^r,>s^. = . (64) 

Note that if only the long-distance El mechanism on the left-hand side of Fig. |6] is con- 
sidered, as was done, for example, in Ref. f25], then current conservation is not satisfied, 
and it appears that some input from short- distance physics is needed in order to define the 
prediction for this observable. 

Since we are considering electric properties, and the spin of the neutron is not affected by 
the photon interaction, it is, once again, convenient to re-express Tg'sfi in terms of the vertex 
function for specific components of the p-wave interaction that appears at the 7r-neutron-core 
vertex: 

^IMr,,. (65) 



2 J 2 ^ 

We note that if we examine the case s' = s then only j = contributes to the sum in 
Eq. (65), and we have 

- r++^. = r__^ = — f 3^ , (66) 
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where ± indicates the spin projection ±1/2. 

For real photons we have k-e = and we choose k-p = 0, with p the incoming momentum 
of the s-wave state. With these restrictions we can write the space-space components of the 
r vertex function as: 



6jiTE + kjPiTM ■ 



(67) 



We can, without loss of generality, choose the photon to be traveling in the z direction, and 
it then follows that = V e- 

Using the definition of B(E1) strength (see, e.g. [3D]) we find that the B(E1) strength for 
this transition is related to the renormalized, irreducible vertex Y ^ by: 



B(E1) 



47r 



r 



++3 



47r V w 



(68) 



where 



1 



(69) 



at leading order. 

Current conservation (64) then provides an alternative way to calculate Te, it tells us 
that: 

ojTjo = kjTE . (70) 
But, for Tjo, the diagram on the right of Fig. [6] need not be considered, and so 



r,o(k) 



Mi(r) 



(71) 



where ui and uq are the leading-order wave functions of the s- and p-wave states, given by 
Eq. (g. As |k| ^ Eq. ^ reduces to: 



Fjo(k) ~ kj / drrui{r) Uo{r) , 



(72) 



an equation in which, of course, the integral is the canonical form of the El matrix element. 
We performed both the momentum-space calculation of F33 = F^, and the calculation 



of F^; via the co-ordinate space integral (72) (with appropriate factors). The divergences 



cancel in the former calculation, as expected, and the same result is obtained from either 



integral. Inserting the result in Eq. (68) we find: 



B(E1) 



72 „2 

^e//e 70 



Stt — ri 

is the leading-order Halo EFT result with 



271 + 7o 
.(70 + 71)^ 



M 



= fQ, ^ 0.366 



(73) 



(74) 



the effective charge. 

No cutoff parameter is needed in order to get a finite result for B(E1): our value is finite 



without regularization, c.f. Ref. [29]. We note that the result (73) is "universal" in the 
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sense that it applies to any El s-to-p-wave transition in a one-neutron halo nucleus. Once 



Ti, 7i, and 70 are known for a given one-neutron halo the prediction (73) is accurate up to 
corrections of order -RcoreZ-Rhaio- 

We now have to deal with the issue that we do not have a value for ri that has been 
obtained solely from an EFT calculation and data. So here we set ri in order to reproduce 
the experimental number ^ obtained in Ref. j3]. This produces: 

r^o = -0.66 fm-\ (75) 

where we do not bother to propagate the error bars from the experiment, since NLO effects 
are presumably a much larger source of uncertainty. This corresponds to Oi = 374 fm^. We 
note, that the value of ri ( [75| ) lies between the two values extracted when we adopt the 
value of ai from Ref. [2J as given and obtained ri using formulae of LO and NLO precision 



in the -RcoreZ-Rhaio expausiou (see Sec. II F). 



Short-distance effects enter B(E1) in these NLO corrections. The B(E1) (1/2^ — )■ 1/2 ) 
transition therefore cannot be predicted at NLO: a counterterm appears at that order. This 



can be seen either from the presence of the operator ~ -^i^i ^ Eq. (|37|) , or from a co- 
ordinate space argument similar to the one made in the previous section for the l/2~ state's 
charge radius. In the case of the El transition we have the co-ordinate space integral |2| 130]: 







^ AA, r dr r e-- ( 1 + ) e"-'^ = 2.p^;°±^ . (76) 

\ lir J V '^i (70 + 71) 



Again, it is finite, and this time it has a short- distance contribution 



^01|5D ~ 2^-^i?core (77) 

which is a fraction ~ -RcoreZ-Rhaio of the total. 

Comparing this calculation with a shell-model treatment of ^^Be, it is clear that one effect 
which is subsumed into the NLO counterterm L^^['^^ is the transition of a neutron from a 
d5/2 to a P3/2 orbital, with that neutron coupled to the 2+ state of ^°Be. This state 
is 3.4 MeV above the ^°Be ground state, so the dynamics associated with it takes place 
at distances -Rcore- Hence in our EFT it can only appear in short-distance operators such 



as that multiplying L^^i'^'' in Eq. (37). The computation of Ref. |1] suggests that such a 



contribution reduces the El matrix element by ~ 30%, which is the anticipated size of an 
NLO effect when the -RcoreZ-Rhaio expansion is employed in the ^^Be system. 

There are other effects of a similar size that will affect B(E1) at NLO. Specifically, there 
are NLO corrections from the wave-function renormalization factors associated with the s- 
wave and p-wave fields. Both tend to increase B(E1) over the leading-order prediction. We 
choose to adjust ri to reproduce B(E1) (1/2"*" — >■ l/2~) already at leading order, and then 
rely on the NLO counterterm to cancel these NLO corrections. That this is a reasonable 
strategy is a testable hypothesis, since the counterterm that appears here at NLO also plays 
a role in Coulomb dissociation of ^^Be. It is to that process that we now turn our attention. 



IV. PHOTODISINTEGRATION OF ^^BE INTO i°BE AND A NEUTRON 

In this section we consider the photodisintegration of the Beryllium-11 nucleus to a 
Beryllium-10 nucleus plus a neutron. In practice this process is measured using Coulomb ex- 



18 



citation of the ^^Be nucleus, with the two reactions connected within the equivalent-photon 
approximation. 

Once again, we first compute photodisintegration into the direct product of a core state 
of definite polarization and a neutron of definite spin. The separation into states of definite 
total angular momentum will be done after that calculation is complete. 

There are two contributions to this process, as depicted in Fig. [7j The first diagram, 
denoted "LO" in the figure, corresponds to the plane-wave impulse approximation contri- 
bution. Evaluation of the relevant Feynman graph yields, for the zeroth-component of the 
four-current: 

,(a) _ eQcgo2mn 



+ (P' - £;k) 

where p' is the relative momentum of the outgoing nc pair. 



(78) 





LO 



NLO 



FIG. 7: The two diagrams representing photodissociation of the ^^Be ground (s-wave) state. As we 
will show below, the first diagram, denoted LO here, is dominant over diagrams involving p-wave 
final-state interactions. 



The second diagram pair of diagrams, (b), includes final-state interactions between the 
neutron and the core. (In fact, these interactions are only present in the J = 1/2 channel, 
a point we will have to deal with below.) The diagram with final-state interactions can be 
written as: 



niR p'^ cot 6i{p') 



with p'^ cot 5i{j)') 



^rip'^ if we do not yet expand for the case ri ~ l/i?core- In 
Eq. ( 79 ) repeated indices are summed over and Lj is the (vector) integral 



ai 



p 



(2vr) 



d-l 



m 



p 



P 



(80) 



This can be converted into a co-ordinate space integral, with the result: 



^ip r 



-Jl 



m 



-uor 



-lor 



Thus, 



M 



p 



12 



2i 



p' -k dr{ip'r - l)e'P''-e-^°' + O 



m 00 
ML ' 



(81) 



(82) 
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Meanwhile, the piece of A^o""* corresponding to a final-state nc p-wave may be rewritten 

as: 

M^^^ = eZeffgo2mRCok -p j dr r^ji{p'r)e-^°'' . (83) 

In the presence of final-state interactions, it is most useful to write the combination A^g"'' -|- 
in the form: 

Aio = eZ(,ffgo2mjiUjk ■ pe^^^^"^ ^ j drr'^ [jiip'r) cos Si{p') + ni{p'r) sin(5i(p')] e~'*'°'' . (84) 

Note that the sign of the spherical Neumann function used here is opposite to that chosen 
in many textbooks. We define: 

, , cosfa;) sinfx) 
ni(x) = ^ + ^^. (85) 



X 



Next we follow Typel and Baur in Refs. [21 [30] and re-express Eq. (84) in terms of 
dimensionless integrals: 



eZe//5fo2m/jw 
where y = p'/To; ^"^^ 



^ effyo n^ f^ . p,^^s,ip') cos{6^{p')) + h{y) sin(5i(p'))] , (86) 

7o 



oo 







My) ^ - I dxx'e-^j,{yx) = -j:^^^^, (87) 



/■°° 1 -\- 3?/2 

h{y)^ dxx\-^n,{yx)= ^ 2,2 2 " (88) 

Jo (1 + y ) y 

This calculation corresponds to the insertion of a plane-wave for the electromagnetic 
field, followed by an expansion of the effects to first order in r. The resulting matrix element 
differs from that of the dipole operator |r|yio(r) by a factor of a/3/ (An) and a l/w. Thus to 
get the physical matrix element of the El strength Mei we multiply A4o by ^y3/{A^^) and 
divide by u. We must also multiply by the field renormalization a/^, for the initial state. 



This also justifies the factors appearing in Eq. (68). 

This produces the El coupling to the physical a field of a plane wave, in the absence of 
neutron spin, J^ei- To get the result for a final state of good total angular momentum we 
need to couple Ai ei to the neutron spinor, and then form states of good total J. This leads 
to two different Mei amplitudes, one for J = 1/2 and one for J = 3/2. 

Only the J = 1/2 amplitude has final-state interaction effects. Thus there we have: 

wi/2) _ eZe//y3 ,,^, ^^^fa,^ 2y3cos(^(i/2)(pO) + (1 + 3y^) sm{6^,/2)ip')) 1 
7o 1/^1 + 1/) V3 

where the irrelevant spin-dependent phase in the Clebsch-Gordan coefficient has been 
dropped. Here we have chosen the photon momentum k to be in the z-direction. The Yio(p') 
that then occurs in A^o is coupled to the final-state neutron spin to produce a 3^(ii)is(p')) 
where s is the spin projection of the initial state. The coupling of Y'io(p') to the neutron 
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spinor yields the Clebsch-Gordon coefficient (|slO|(^l)^s) = The projection onto the 

J = 1/2 continuum nc state then removes this 3^. Note that physically the final-state spin 
projection of the core-neutron state is equal to that of the initial state because the of 
the neutron is not changed by the electric transition. Meanwhile, the angular-momentum 
projection of the nuclear state is unaffected by the photon, since it is taken to be in the 
z-direction. Finally, the final-state strong interactions conserve the sum rris + mi, and so m'j, 
the angular-momentum projection of the final state is equal to s, the spin projection of the 
neutron in the initial state. 

For the spin-3/2 part final-state interactions are natural, and so are suppressed by three 
powers of i?core/-Rhaio- Thus in this piece of the amphtude we neglect diagram (b), which is 



equivalent to taking 5i — > in Eq. (86). Alternatively, we can project out the El, J = 3/2, 
piece of the result from Eq. (78) for and multiply by y/Z^ to find: 



= ^oeZ,,,y3^-;,^^, (90) 

where an irrelevant spin-dependent phase again has been dropped. In this case the coupling 
of Yio{p') to the neutron spinor yields a Clebsch-Gordan coefficient (^slO|(^l)|s) = 
Finally we convert to the physical observable ^^j^^^ , by: 

dB(El)=(|ATri=+l>fr'lOjgi. (91) 



Inserting Eqs. (89) and (|90|) in Eq. (91) produces: 



rfB(El)_ 2 mn,2f W cos{6^i/2){p')) + {l! + S^op") sm{5^,/2){p'))V , 8p'' 



dE ^^^27r2 0^ p'3(p'2 + 72)4 (p'^ + li)^ 

(92) 

where E = p'"^/ {2m r) is the kinetic energy of the neutron and ^°Be in the center-of-mass 
frame. This expression can be rendered more straightforward to work with if we use the 
identity cot^ 6 + 1 = sin~^ 6 to write things in terms of p'^ cot 5(i/2)(p')- 

c^B(El) _ 2 mn ^ [ p'W cot (^(1/2) (pQ) + 7^ + ^lop"? 8p'' \ 

dE ^^^27r2 ° V [p"+P"cot2(5(i/2)(p'))](p" + 7o)' (p" + 7o)V ' 

which makes it clear that the first term (the contribution of the J = 1/2 channels) does not 
diverge as — )• 0. 



The expression (93) is true to all orders in final-state interactions, and for any value of 
the s-wave asymptotic normalization Aq. However, as explained above, in the ^^Be system, 
final-state interactions — even in the J = 1/2 channel — are weak, and can be considered 
sub-leading. Also, Aq = \/2'Jq at leading order, but then receives higher-order corrections in 
the case of a non-zero tq. Thus we now provide order-by-order expressions for (iB(El)/dE 
in the Rcore/ -Rhaio expausiou. 

First, at LO we have 5(i/2)(p') = 0, and Aq = y/2jQ, so: 

^B(E1) ^"_ , , 3mn Sjof' 
dE ^ ^^^27r2 (p'2 + ^2)4- y^'V 
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The NLO correction comes from two sources. The first is the shift of Aq to larger values 
due to a To > 0, which tends to increase the B(E1) strength. Second, final-state interactions 



make 5(i/2)(p') 7^ 0. Employing Eq. (32) for 5(1/2) (p') in Eq. (92) and keeping only terms up 



to first order in that phase shift, the total NLO result becomes: 



dB(El) ^^° _ 2 3m« 8^op'-' ( Al 2^o ^l + ^p'^ \ 
dE 27r2 (p'2 + ^2)4 1^270 ^ 3ri p'^ + 7? ) 

SrriR 8jop'^ f 270 jl + 3^\ 



272 3mij 



where in the second line we have used ( 18 ) and kept only the term of 0(ro7o). Accurate mea- 
surements of the Coulomb dissociation spectrum therefore provide information on the s-wave 
^°Be effective range, if the p-wave effective range is already fixed from another observable. 



This sort of expansion of Eq. (93) can be carried to arbitrarily high orders in i?core/-Rhaio- 
However, it must be remembered that the counterterm that affects B(El)(l/2+ — t- 1/2^) 
will appear in dB{El) / dE at 0{Rl^j.^/Rlg^i^), and so the expansion of FSI and asymptotic- 
normalization effects in powers of the EFT small parameter does not capture all of the 
physics of the B(E1) transition. 

We do, however, have the advantage that the same counterterm enters both the bound- 
to-bound and bound-to-continuum transition, so the number of free parameters needed to 
describe B(E1) strength is limited. Up to NLO accuracy for the bound-to-bound, and 
NNLO for the bound-to-continuum, there are five numbers: 70 and 71 (which are known 
from separation energies), and ri, Aq, and the counterterm, which must be fitted to the 
available data on B(E1) (1/2+ -> 1/2") and dB(El)/dE. 

The quality of the Coulomb excitation data is such that realistically we can hope only 
to extract one parameter from it, and so an NNLO analysis is not feasible at this time. In 
the next section we present our numerical results from the NLO analysis of the Coulomb 
excitation data. At NLO we need only 70, 71, ri, and Aq (or, equivalently, tq). 



V. RESULTS FOR COULOMB EXCITATION OF THE "BE NUCLEUS 



The data set we use is that of Ref. |5j. The predictions obtained in the previous section 
must be folded with the neutron detector resolution as described in that paper. For the 
number of El photons as a function of the photon energy u we take 



Nei{co) 



iKo (0 K, (0 - [Kl (0 - Kl (0] 



(97) 



where ^ = uh/ ("fv), with b the impact parameter, and 7 = (1 — t'^/c^)"^/^ ^^j-^g Lorentz factor 
of the beam. We adopt the value b = 10.38 fm and obtain 7 from the kinetic energy of the 
^^Be beam of 520 MeV per nucleon, both based on values quoted in Ref. 



The leading-order result ( 94 ) can then be folded with detector response and the spectrum 
of El photons. This produces the red-dashed curve in Fig. |8| At next-to-leading order we 
take the value of ri fixed as per Eq. (75). We then have one free parameter, the value of Aq 
at NLO (or, equivalently tq) [32]. After folding we find a reasonable fit for Aq/A^'^ = 1.3, 
which corresponds to tq = 2.7 fm. The green dotted line shows the result of Typel and Baur 
of Ref. [30]. That result was obtained with all integrals regulated at a scale of i? = 2.78 
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FIG. 8: Differential B(E1) strength for Coufomb dissociation of Beryllium- 11 into a neutron and 
a ^'^Be nucleus, plotted versus the excess energy of the detected neutron E*, in MeV. The data 
are from Ref. 0. The theory curves have been folded with detector resolution. The red dashed 
line is the leading-order Halo EFT prediction, which does not include any final-state interactions. 
Final-state interactions, with the effective range taking on a value fixed from the bound-to-bound 
El transition strength, are included in the NLO result, which is shown in blue. The result of 
Ref. [30] is the green dotted line, which essentially matches the solid blue line. 

fm. This corresponds to specific assumptions about all the counterterms that appear in the 
theory. 

Note that Aq/Aq'~' = 1.3 increases the charge radius of the ^^Be ground state to 

{rl)[ll = 2.42 fm, (98) 

which is certainly in agreement with the atomic-physics number ([s]) within the expected size 
of NNLO corrections. This must be taken with a grain of salt, though, since the change in 
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the A = 10 to A = 11 radius difference from LO to NLO is ~ 30%. But, in contrast to 
observables involving p-wave binding, the radius of the ^^Be ground state does not receive 
any corrections from short-distance physics until N^LO. The remaining difference between 



our NLO number (98) and the experimental number is certainly consistent with the presence 



of the short-distance operator ~ L^(jQ &t N'^LO in the EFT expansion for the radius. 

The neutron radius of the ground state of ^^Be is similarly increased by the presence 
of an ro > 0. The value of Aq extracted from the Coulomb dissociation data produces 
(r^)^/^ = 5.6 fm within the context of our NLO calculation, a 30% shift from the LO result. 
As with {r\) , the next contribution is from purely short-distance physics, and assuming that 
this gives a contribution of order -Rcorc to (rn)^^"^, we have the Halo EFT determination: 

(r2)V2 = 5 ± 0.6) fm, (99) 
where the error bar does not account for the statistical uncertainty in the extraction of Aq. 



Eq. (99) is a model-independent result for the neutron radius of ^^Be, obtained from the 



Coulomb dissociation data of Ref . [3] . 



VI. CORRELATIONS BETWEEN OBSERVABLES IN THE ^ BE SYSTEM 



Effective field theories in general, and Halo EFT in particular, provide model-independent 
correlations between different observables. In the previous sections, we have expressed the 
electromagnetic properties of the ^^Be system through the effective range parameters for 
n-^°Be scattering, 70, 71, tq, and ri. Our expressions thus can be interpreted as correla- 
tions between scattering observables and electromagnetic properties. However, depending 
on the experimental information available, it could be useful to look at correlations between 
electromagnetic observables. 

As an example, we consider the correlation between the B(E1) strength and the radius 
of the p-wave state in ^^Be at LO. Using Eqs. (58, 73) we obtain 



B(E1) 



15n 



2\W 



X 



l + 2x 

(1 + X)2 



(100) 



where x = ^^/Bi/Bq is the ratio of the neutron separation energies for the p-wave and s-wave 
states. The B(E1) strength is thus proportional to the mean-square radius of the p-wave 
state. In the limit of vanishing neutron separation energy for the p-wave state, the B(E1) 
strength vanishes linearly with x. Eq. (100) can also be used to obtain the electric radius 
of the p-wave state {r'^)^'^^ directly from the measured value of B(E1) and the neutron 
separation energies Bi and Bq. This gives a radius of the ^^Be p-wave state, relative to the 
lOgg ground state: 

(r|)W = 0.35.. .0.39 fm^ (101) 

depending on whether the value for B(E1) from Eq. ^ or Eq. ^ is used. 

Analogously, the strength for the El transition to the continuum can be related to the 
radius of the s-wave state. However, this observable is also affected by the magnitude of 
p-wave final-state interactions. At next-to-leading order we have: 



dB(El) 
dE 



NLO 



12 



y 



n^Bo (1 + 2/2)4 



e^Q^A(r|)('^) -7rB(El) 



{l+x)\l + 3y^) 
{y"^ + x^){l + 2xf 



(102) 



24 



where y = p'/To- Similar correlations can be derived for other observables and/or at higher 
orders. They have proven useful in the analysis of universal properties in ultracold atoms 
and also show promise for halo nuclei [33\ . 

VII. CONCLUSION 

This discussion of electromagnetic observables in the Beryllium- 11 system already dis- 
plays the significant recent experimental activity that has been focused on this system. 
Coulomb excitation has been used at a variety of facilities to probe El transitions, and 
atomic-physics experiments have made great strides through advances in trapping technol- 
ogy. This means that the time is ripe for a detailed analysis of electromagnetic properties 
of halo systems. Here we have shown how EFT can provide such an analysis. 

The Halo EFT we employed is complementary to ab initio methods [T014T2] . which can 
struggle to describe El transitions and radii in these extended nuclei because of the widely 
varying core and halo scales that are present in the problem. In Halo EFT this wide 
separation of scales is the basis for the calculation. Input that summarizes the physics 
at scale -Rcore can be taken from either simulation or experiment, and the EFT is then used 
to predict the outcome of experiments that probe dynamics at the halo scale. 

The results derived here in EFT mirror elegant analytic approaches to halo nuclei (see, 
e.g. Refs. [ZlEniEO]). But, in contrast to those works, there is no regulator dependence in our 
result for the El strength. This is a consequence of current conservation in our formalism. 
Moreover, the EFT delineates the order at which any observable receives a contribution 
from physics at scale -Rcore — a contribution which cannot be calculated using the asymptotic 
wave functions ([6]). The EFT also explicitly shows whether — and if so, where — these short- 
distance effects appear in other observables. In this way we can systematically assess the 
impact of physics at scale -Rcore on low-energy electromagnetic observables in halo nuclei, 
and so go beyond the calculations of Refs. EOj. 

We determined the magnitude of the s-wave n-^°Be effective range, tq, by examining 
experimental results for the low-energy El strength function in breakup to the ^''Be-neutron 
channel: dB{El)/dE. We find a reasonable fit for ro = 2.7 fm. We were also able to extract 
a value for the p-wave effective range ri from the bound-to-bound B(E1) strength, obtaining 
ri = —0.66 fm~^ up to the 40% corrections which seem typical in this EFT. This gives a 
p-wave scattering volume ai = (374 ± 150) fm'^. This overlaps the range quoted in Ref. [30] , 
but has a lower central value. We would caution against any strong conclusion regarding 
values of ai and ri until the counterterm that parameterizes the short-distance piece of the 
El strength in this system has been determined by further measurements. 

There are no spectroscopic factors in our approach. Nevertheless, a comparison of the 
values we have obtained for the physical observables ro, ri and ai (or, equivalently, Aq 
and Ai) with those found within many-body models and ab initio calculations would be 
interesting. While the interpretation of these quantities in terms of the many-body dynamics 
of the A = 11 system is not one of the goals of this work, such an interpretation could inform 
which observables are interesting for future investigation within the EFT, and, in particular, 
which ones will most reward study at higher EFT orders. 

In this vein, we note that a higher-order calculation of the El observables considered here 
would necessitate determination of the El counterterm that occurs at NLO in the bound- 
to-bound B(E1) strength. Better data on Coulomb dissociation of ^^Be at low energies 
would be very helpful in this regard, since the same counterterm occurs in the continuum 
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El strength at next-to-next-to-leading order. With this one additional parameter we should 
be able to perform calculations to 5% or better for dB(El)/dE. It is, though, important 
to note that the EFT used here breaks down at energies of order a few MeV, and so the 
window in which the theory will give such an accurate description of data is quite limited. 

But, already at the order to which we have worked here, the extraction of parameters for 
n-^°Be scattering facilitates predictions for the electric and neutron radii of ^^Be. We find 
agreement at the level expected of our NLO calculation between our result for ('"l;)iiBg and 
the measurement of Ref. [8]. We also predict that the excited state of ^^Be has an electric 
radius of 2.43 ± 0.1 fm. Meanwhile, the neutron radius of the ^^Be ground state is found to 
be (r^)^/^ = 5.6 ± 0.6 fm. We note that EFT provides a way to extract this quantity from 
the Coulomb dissociation data which is independent of any assumptions about the details 
of the physics that occurs at distances ~ -Rhaio- 

Lastly, we emphasize that the Halo EFT can be used to derive correlations driven by 
the separation of scales in the system. These correlations may not be obvious in ab initio 
approaches where this separation is not explicit in the parameters of the theory. EFT 
has been successful at revealing these patterns in few-nucleon systems, e.g. the correlation 
between the triton binding energy and the spin-doublet neutron-deuteron scattering length 
(the "Phillips line" ^Ij), and the correlation between the triton and alpha particle binding 
energies (the "Tjon line" ^^'^ understand that those correlations are driven by 

the hierarchy between NN scattering lengths and the range of the NN force: l/a^, 1/at <C 
m^r, and so they can both be derived within an appropriate EFT EZ]- Here we have 
presented analogous correlations for electromagnetic observables in the ^^Be system. Similar 
relations have proven of significant interest in the field of ultracold atoms [33] . Future tests 
of correlations such as Eqs. (101) and (102) can reveal the extent to which the physics of 
^^Be, and that of other halo nuclei, is driven by "universality". 

The application of the Halo EFT to other one-neutron halos, and to two-neutron halos 
such as ^^Li, are obvious next steps. With the inclusion of Coulomb interactions, proton 
halos become accessible as well. 
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